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A FINITE REFLECTION FORMULA FOR A POLYNOMIAL 
APPROXIMATION TO THE RIEMANN ZETA FUNCTION 


STEPHEN CROWLEY 


ABSTRACT. The Riemann zeta function can be written as the Mellin transform of the unit interval 
map w(x) = [271] (x [271] + 2-1) multiplied by sH. A finite-sum approximation to ¢ (s) 
denoted by Çw (N;s) which has real roots at s = —1 and s = 0 is examined and an associated 
function x (N; s) is found which solves the reflection formula Gw (N;1—s) = x (N; s) Cw (N; s). 
A closed-form expression for the integral of Cw (N; s) over the interval s = —1...0 is given. The 
function x (N; s) is singular at s = 0 and the residue at this point changes sign from negative to 
positive between the values of N = 176 and N = 177. Some rather elegant graphs of Cw (N; s) 
and the reflection functions x (N; s) are also provided. The values w (N;1 — n) for integer values 
of n are found to be related to the Bernoulli numbers. 
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1. THE RIEMANN ZETA FUNCTION AS THE MELLIN TRANSFORM OF A UNIT INTERVAL MAP 


The Riemann zeta function can be written as the Mellin transform of the unit interval map 
w(x) = |] (x |a-1| + 2-1) multiplied by s$. [3][2] 
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FIGURE 1. The Harmonic Sawtooth map 


1.1. The Truncated Zeta Function. The substition co — N is made in the infinite sum appear- 
ing the expression for („(s) to get a finite polynomial approximation 
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with equality in the limit except at the negative integers 
(3) limysooGu(N;s) = Cl) sg N* 


and where Y (x,n) = 34, V (x) is the polygamma function and Y (x) = In (T (z)) is the digamma 
function. The functions Çw (N; s) have real zeros at s = —1 and s = 0, that is 


(4) dim, Cu (N; s) = lim Cu (N;s) =0 
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One possible idea is that the functions („(N;s) can be orthonormalized over the interval s = 
—1...0 via the Gram-Schmidt process[4] and that the result might possibly shed some light on the 
zeroes of ¢(s). Let the logarithmic integral be defined 
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where y = 0.57721... is Euler’s constant, then the normalization factors are given by the integral 
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FIGURE 2. [Cw(N;s):s=-—1...0,N =1...25} 


The following table lists the values of ¢,, (N; 1 — n) for n = 2...12. 
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4 CROWLEY 
1.1.1. Integrating Over the Critical Strip. There is a formula similiar to (6) which gives the integral 
of Cw (N; s) over the critical strip 0 < Re (s) < 1. 

J Cw (N;c+is)de =1+ Wet (Ei (isla (N + 1) — ln (N + 1)) - Ei (isn (N + 1))) 


N-1  n(n+1)7*8 
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(7) 
where Ei; (t) is the exponential integral defined by 
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The contribution from the Ei term vanishes as s — oo, that is 


; N LE de. 5 
(9) lim al (Ei (ts In (N +1) — ln (N + 1)) — Ei (is In (N +1))) =0 


1.1.2. The Reflection Formula. There is a reflection equation for the finite-sum approximation 
Cw(N; s) which is similiar to the well-known formula ¢ (1 — s) = x (s) € (s) with 


E TS 
x (5) = 2(2m)~* cos (=) T (5). 
The solution to 
(10) Cw (N;1— 8) = x (N; 5) Cw (N; $) 
is given by the expression 


x(N;s) = Sw(N;1=s) 
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which satisfies 

(12) x(N;1-— 8) = x (N; s)! 

The functions x (N; s), indexed by N, have singularities at s = 0. Let 
a (N) = Y, n (In (n +1) — In (n)) 

(13) b(N) = yi In(n)n? —In(n+1)n? —In(n) 


n=l n(n+1) 


c(N) =1 5, m (In (n +1)? — In (ny) 








then the residue at the singular point s = 0 is given by the expression 
(14) 
Res(x(N;s)) = —Res(x(N;8)~") 
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which has the limit 
(15) lim Res(x(N;s)) =1 


No s=0 


A FINITE REFLECTION FORMULA 5 


We also have the residue of the reciprocal at s = 2 
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which vanishes as N tends to infinity 
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As can be seen in the figures below, the residue at s = 0 changes sign from negative to positive 
between the values of N = 176 and N = 177. 
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FIGURE 3. {Res(x(N; s)):N= 1...250} 


For any positive integer N, we have the limits 


lims>0 x (N; s) = 00 

lims_50 SEX (N;s) =œ 

lim, 1 x (N;s) =1 

18 s> Fj 
(18) lim X (N; s) =0 
lims_,2 X (N; s) =0 
lim, 1 Lx (N;s) =0 


The line Re (s) = 3 has a constant modulus 


1 
(19) x (mig +s) 1 
There is also the complex conjugate symmetry 


(20) X(N; a + ty) = x (N;x — iy) 
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FIGURE 4. [Res(x (N; 5)" :N=1... 250} 


If s =n € N* is a positive integer then x (N;n) can be written as 
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where (,, ,) is of course a binomial. e Bernoulli numbers|1| make an appearance since 
here (71) is of b l. The B ll bers[1] mak pp 


(22) X(N32n) Go(N32n) = Ban (N +1)? 2D 4... 
The denominator of x (N;n) has the limits 
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Another interesting formula gives the limit at s = 1 of the quotient of successive functions 
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FIGURE 
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Let 
(25) v (s) = x (09; 5) = ~ 


Then the residue at the even negative integers is 


¢(1=n) 
Dr: Sy DO ES neven 
26 Res (v(s)) = #¢(s)|s=—n 
(26) Res (vis))=4 ¿ BR 
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